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Abstract. In 1952 Lee and Yang proposed the program of analyzing phase 
transitions in terms of zeros of partition functions. Linear operators preserving 
non-vanishing properties are essential in this program and various contexts in 
complex analysis, probability theory, combinatorics, and matrix theory. We 
characterize all linear operators on finite or infinite-dimensional spaces of mul- 
tivariate polynomials preserving the property of being non-vanishing whenever 
the variables are in prescribed open circular domains. In particular, this solves 
the higher dimensional counterpart of a long-standing classification problem 
originating from classical works of Hermite, Laguerre, Hurwitz and Polya- 
Schur on univariate polynomials with such properties. 
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Introduction 

Zero loci of multivariate polynomials or transcendental entire functions and their 
dynamics under linear transformations are central topics in geometric function the- 
ory which in recent years has found applications in statistical mechanics, combina- 
torics, probability theory and matrix theory |6l HI [TOl HH [El El 113 • For instance, 
as shown in [6l [131 132] , the Lee- Yang program on phase transitions in equilibrium 
statistical mechanics (31] [59] is intimately related to the problems of describing 
linear operators on polynomial spaces that preserve the property of being non- 
vanishing when the variables are in prescribed subsets of C". In this paper we 
fully solve these problems for (products of) open circular domains. In particular, 
this accomplishes a long-standing classification program originating from classical 
works of Hermite, Laguerre, Hurwitz and Polya-Schur on univariate polynomials 
with such properties [H [H [H SIl 111113] ■ 

Given an integer n > 1 and H C we say that / £ C-lzi, . . . , Zn] is fl-stable 
if f{zi,...,Zn) =/= whenever (zi,...,2:„) S f2. A K-linear operator T : V 
K[zi, . . . ,Zn], where IK = K or C and ^ is a subspacc of K[zi, . . . , z„], is said to 
preserve ^-stability if for any fi-stable polynomial f (z V the polynomial T(f) is 
either fi-stable or T{f) = 0. For k = . . . , k„) e N" let Kk[2;i, ■ • ■ , ^n] = {/ G 
K[zi, . . . , Zn] : deg^. (f) < Ki,l < i < n}, where deg^. (/) is the degree of / in z^. If 

C C and n = then il-stable polynomials will also be referred to as ^f-stable. 
The fundamental classification problems that we address below are as follows. 

Problem 1. Characterize all linear operators T : K„[zi, . . . , z„] ]K[zi, . . . , z„] 
that preserve Cl-stability for a given set C C" and k G N". 

Problem 2. Characterize all linear operators T : Klzi, . . . , Zn] — > K[2i,...,x„] 
that preserve fl-stability, where Q, is a prescribed subset o/C". 

For n — 1 these problems were stated in this general form in [TH [T7] , see also 
[25] and [131 pp. 182-183]. Their natural multivariate analogs given above were for- 
mulated explicitly in [9l [10] , thereby encompassing essentially all similar questions 
or variations thereof that have been much studied for more than a century. Note 
that for n = 1, K = K and = {z e C : Im(z) > 0} Problems [l]-[2| amount to 
classifying linear operators that preserve the set of real polynomials with all real 
zeros. This question has a long and distinguished history going back to Laguerre 
and Polya-Schur, see [10l[T4] and the references therein. In particular, in ^4T] Polya 
and Schur characterized all diagonal operators with this property (so-called mul- 
tiplier sequences), which paved the way for numerous subsequent investigations 
[a[9l[l5l[l6l[T7l[ll[26l[32l[35l[40l[43l[5^ However, it was not until very recently 
that full solutions to this question - and, more generally, to Problems [T]-[5| for n ~ 1 
and any open circular domain f2 - were obtained in [7] . 

In physics literature ([51], [53l §5]) one sometimes distinguishes between "soft" 
and "hard" theorems asserting the non-vanishing of partition functions in certain 
regions. This dichotomy stems from "soft-core" vs. "hard-core" pair interactions in 
lattice-gas models and does not refer to the level of difficulty in proving such theo- 
rems but rather to the fact that in some sense "soft" theorems are constraint-free 
while "hard" theorems involve constraints of various kinds, such as the maximum 
degree of a graph. By analogy with this terminology, one may say that results per- 
taining to Problem [1] are "hard" or "algebraic" (bounded degree) while those for 
Problem [2| are "soft" or "transcendental" (unbounded degree). For n > 1 several 
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partial results of the latter type were established and appplied to various questions 
in e.g. combinatorics and statistical mechanics in the past two decades. In j33j Lieb 
and Sokal proved a general result [331 Proposition 2.2] pertaining to Problem [2] and 
more results of the same kind were recently obtained in ^13l I57j . Another contri- 
bution to this subject is Hinkkanen's Schur-Hadamard composition theorem [23j . 
Progress on Problem [2] was recently made in 9J where a complete solution was 
given for finite order partial differential operators when Q = ^E*" and ^' is an open 
half-plane. Further contributions to this problem have been reported in ,19J. 

The main results of this paper provide complete solutions to Problems [T]-[5] when 

— Qi X ■ ■ ■ X Qn and ^i, 1 < i < n, are arbitrary open circular domains in 
C. We also characterize all linear transformations preserving the class of Lee- Yang 
polynomials defined with respect to any such set Q. In particular, the classification 
theorems of |7J may now be viewed as special cases of the theorems below from 
which they follow by setting n = 1. To achieve this we need several new ideas and 
results. For instance, we define the polarization of a linear operator to reduce the 
sufficiency part for arbitrary degrees in the algebraic characterization to the case of 
linear transformations acting on polynomials of degree at most one in each variable. 
For the necessity part in the transcendental characterization we generalize Szasz' 
inequalities to several variables. These are bounds on the coefficients of a stable 
polynomial that only depend on its first few non-zero coefficients. 

The solutions to Problems [THl] may be summarized as follows: (essentially all) 
linear operators preserving n-variate stable/Lee- Yang polynomials are induced by 
2n-variate stable/Lee- Yang polynomials via appropriately defined symbol maps. 

In a sequel f6| to this paper we build on our classification theorems to develop a 
self-contained theory of multivariate stable polynomials. We therefore take care to 
rely on as few auxiliary results as possible in the process. Combined with the present 
work, the theory and applications in ^ also yield a natural framework for dealing in 
a uniform manner with Lee- Yang type problems in statistical mechanics, combina- 
torics, and geometric function theory in one or several variables, thus contributing 
to Ruelle's quest |46] for an appropriate mathematical context encompassing the 
celebrated Lee- Yang theorem [31( (locating the zeros of the partition function of 
the ferromagnetic Ising model on the imaginary axis in the complex fugacity plane) 
and its many modern relatives [ll[3ll[5l[2a[2a[33[3a[MlllIlllHll5l]. This also 
illustrates Hinkkanen's observations about zeros of multivariate polynomials and 
their "so far unnoticed connections to various other concepts in mathematics" [2 3) , 
which are further substantiated by recent applications to probability theory, matrix 
theory and combinatorics [8llITl[12iri3, 34, 44, 54, 57J. 

1. Operator Symbols and fi-STASLE Polynomials 

The general notion of f2-stability defined in the introduction extends classical 
univariate concepts such as Hurwitz or continuous-time stable polynomials and 
Schur or discrete-time stable polynomials (see, e.g., [351133]). These correspond to 
n = I, n = {z £ C : Re{z) > 0} and n = 1, = {z e C : 1^1 > 1}, respectively. 

Let now n be an arbitrary positive integer. For 6 £ [0, 2tt) set 

He = {z e C : Im(e^^z) > 0}. 
We find it convenient to work with the upper half-plane 

H -.^Mo^izeC: lm{z) > 0}, 
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therefore we will follow Levin's terminology [32] and refer to i/-stable multivariate 
polynomials simply as stable polynomials (cf. [3 El [9l [11] ) . A stable polynomial 
with all real coefficients is called real stable. Clearly, a univariate real polynomial is 
stable if and only if all its zeros are real. We denote the sets of stable, respectively 
real stable polynomials in n variables by 7i„ (C) and 7i„ (R) , respectively. 

If = M:^ = {z & C : Re(z) > 0} then il-stable polynomials are called weakly 
Hurwitz stable. In I13j these are termed polynomials with the half-plane property. 

1.1. Algebraic and Transcendental Symbols. The solutions to Problems [THD 
for circular domains make use of appropriately defined operator symbols that we 
proceed to describe. For simplicity, in this section we only state the main theorems 
for linear operators preserving (_ff-)stability. 

Recall that Kk[2:i, . . . , 2;„], where K = R or C and k ~ . . . , k„) e N", is the 
K-space of polynomials in variables zi, . . . , z„ of degree at most Ki in Zi, 1 < i < n. 
If 7 = (71, . . . , 7„) e N" we define the algebraic symbol of a linear operator 

T : K„[zi, . . . , z„] ^ lK^[zi, . . . , Zn] 

to be the polynomial Gt{z, w) G K^©k[-Zi, • ■ • , 2„, wi, . . . , w„] given by 

Gt{z,w)^T[{z + wY] - {'^)t{z'^)w^~'^, 

where jOk = (71, . . . , 7„, ki, . . . , k„) while < denotes the standard product partial 
order on N" and we employ the usual multi- index notations z = (zi, . . . , z„), w = 
(wi, . . . ,w„), z° = nr=i = (ziwi, . . . ,z„u;„) and 

\ctj lo otherwise. 

Our first theorem provides an algebraic characterization of stability preservers 
on finite-dimensional complex polynomial spaces, which solves Problem[T]for K = C 
and fl = H. 

Theorem 1.1. Let k E N" and T : Ck[zi, . . . , z„] — + C[zi,...,z„] be a linear 
operator. Then T preserves stability if and only if either 

(a) T has range of dimension at most one and is of the form 

T{f) = a(/)P, 

where a is a linear functional on Ck[zi, . . . , z„] and P is a stable polynomial, 
or 

(b) GT{z,w)En2n{'C). 

Remark 1.1. Theorem 11.11 mav also be stated in terms of other (algebraic) symbols 
of T that are equivalent to the one defined above, for instance the polynomial 

a<K 

Note that T[(l - zwY] G 7i2«(C) ^ Gt{z,w) G 7i:2„(C), see, e.g., Lemma[r7|(3). 

Definition 1.1. Two polynomials /, g G M[zi, . . . , z„] are in proper position, de- 
noted f <^ g, it g + if e 7i„(C). 
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By the Hermitc-Biehler theorem (see, e.g., [43l Theorem 6.3.4]) in the univariate 
case the relation f ^ g is equivalent to saying that / and g are real-rooted (or 
identically zero), their zeros interlace and their Wronskian W[f,g] := f'g — fg' is 
non-positive on the whole of M. We study further properties of this relation in i )1.2l 

The algebraic characterization of real stability preservers on finite-dimensional 
polynomial spaces - which solves Problem [T] for K = R and — H - \s as follows. 

Theorem 1.2. Let k G N" and T : R^i^i, • ■ ■ , 2„] R[zi, . . . , 2;„] he a linear 
operator. Then T preserves real stability if and only if either 

(a) T has range of dimension no greater than two and is of the form 

T{f)^a{f)P + l3{f)Q, 

where a, (3 : Rk[zi, • . • , ^^n] ^ R o,re linear functionals and P,Q are real 
stable polynomials such that P <^ Q, or 

(b) Gr(z,w) eH2„(R), or 

(c) Gt(^,-H e7i2„(M). 

If T : K[zi, . . . , Zn] ~* ]K[zi, . . . , z„], where K = M or C, is a linear operator 
we define its transcendental symbol, Gt{z,w), to be the formal power series in 
Wi, . . . , Wji with polynomial coefficients in ]K[zi, . . . , z„] given by 

GTiz,w)= ^ (-l)"T(z")^, 
where a\ — ail ■ • • a„!. By abuse of notation we write Gt{z, w) — r[e~^ "'], where 

Z ■ W = ZiWi -I- ... -I- ZnWn- 

Let us define the complex Laguerre-Polya class 7i„(C) as the class of entire 
functions in n variables that are limits, uniformly on compact sets, of polynomials 
in 7i„(C), see, e.g., \3T, Chap. IX]. The usual Laguerre-Polya class 7i„(R) consists 
of all functions in 7i„(C) with real coefficients. 

Our next theorem provides a transcendental characterization of stability pre- 
servers on infinite-dimensional complex polynomial spaces, thus solving Problem [2] 
for K = C and H. 

Theorem 1.3. Let T : C[zi, . . . , z„] — s- C[zi, . . . , z„] be a linear operator. Then T 
preserves stability if and only if either 

(a) T has range of dimension at most one and is of the form 

T{f) = a(/)P, 

where a is a linear functional on C[zi, . . . , z„] and P is a stable polynomial, 
or 

(b) Gt(z, w;) eH2„(C). 

Remark 1.2. From Theorem 11.31 one can easily deduce a characterization of linear 
operators preserving fJ-stability for any open half-plane 17. For instance, the analog 
of Theorem 11.31 (b) for the open right half-plane is that the transcendental 
symbol ofT with respect to Hz., i.e., the formal power series 

T[e--] := ^ r(z")^, 

defines an entire function which is the limit, uniformly on compact sets, of weakly 
Hurwitz stable polynomials. 
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The analog of Theorem 1 1.31 for real stability presevers - which solves Problem [2] 
for K = R and fl ^ H ~ is a.s follows. 

Theorem 1.4. Let T : R[zi, . . . , 2„] — > . . . , z„] be a linear operator. Then T 
preserves real stability if and only if either 

(a) T has range of dimension no greater than two and is of the form 

T{f)=a{f)P + (3{f)Q, 
where a, /3 : M[zi, . . . , z„] ^ M are linear functionals and P, Q are real stable 
polynomials such that P Q, or 

(b) Gt(2,w) e7?2^(M), or 

(c) GT(2,-t«) GH2„(M). 

1.2. Fundamental Properties of Stable Polynomials. In this section we first 
review some basic facts and then prove several results on (real) stable polynomials 
needed later on. We start with an immediate consequence of the definitions (cf. [9l 
Lemma 2.1]). 

Lemma 1.5. Let f G K[2;i, . . . , Zn], where IK = R or C. Then f G 7i„(K) if and 
only if f{Xt + a) e HiiK) for all A e M.^ and a G R". 

The following theorem is a multivariate version of Hurwitz' theorem on the 
"continuity of zeros", see, e.g., [131 Footnote 3, p. 96] for a proof. 

Theorem 1.6 (Hurwitz' theorem). Let D be a domain (open connected set) in C" 
and suppose {fk}kLi ^ sequence of non-vanishing analytic functions on D that 
converge to f uniformly on compact subsets of D. Then f is either non-vanishing 
on D or else identically zero. 

We next list some of the closure properties for (real) stable polynomials. The 
first property in Lemma 11.71 below is deduced by applying Theorem 11.61 with D — 
{z e C : Im(z) > 0}" and /fc(zi, . . . , z„) = f{zi, Zi-i, ^ + Zi/k, z^+i, . . . , Zn), 
fc e N, while the remaining properties are easy consequences of the definitions and 
Lemma 11.51 

Lemma 1.7. Let K = R or C and f E 7Y„(K) be of degree dj in Zj , 1 < j < n. 
Then for any 1 < i < n one has: 

(1) f{zi,--- ) e Hn-i{K) U {0} for fi e R; 

(2) /(zi,...,2;j_i,Azj,Zi+i,...,z„) e 7Y„(IK) for X > 0; 

(3) zf'f{zi, . . . , -z^"\ Zi+i, . . . , z„) e Hni^); 

(4) /(zi, . . . ,Zi_i,Zj, Zi+i, . . . ,z„) e 7i„_i(IK) for I < j ^ i <n. 

We will now establish a series of results involving the notion of proper position 
introduced in Definition ll.il compare with [Q, §2], |13i, Theorem 3.2]. In particular, 
the next lemma provides a multivariate analog of the classical Hermite-Biehler 
theorem showing that proper position is a natural higher dimensional counterpart 
of interlacing. 

Lemma 1.8. Let fjgE R[zi, . . . , z„] \ {0}, set h — f -\- ig and suppose that f and 
g are not constant multiples of each other. The following are equivalent: 

(1) h e TiniC), that is, g < f; 

(2) \h{z)\ > \h{z)\ for all z = (zi, . . . ,z„) G C" with lm{zj) > 0, I < j < n, 
where z = (zi, . . . , z„); 
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(3) f + Zn+igeHn+ii^); 

(4) /,gG H„(R) and 



whenever z = (zi, . . . , Zn) € C" with lm{zj) > 0, I < j < n; 
(5) g{Xt + a)<^ f {Xt + a) for all A e R'^ and a e M". 

Proof. The equivalences (1) <^ (5) and (3) <^ (4) are immediate from the definitions 
and Lemma fL5l It is also clear that (3) =i> (1) and (2) ^ (1). We proceed to prove 
the remaining implications, namely (f) => (3) and (f) (2). 

Suppose that (1) holds, i.e., h S Hn(C), and w ~ a + i(3 £ {( e C : Ini(C) > 
0}" with a, /3 e R". The polynomial p{t) = h{a + t(3) is stable so we may write it 



p{t) = Cll{t-Q), 

where C e Candlm((^j) < 0, 1 < j < d. Note that for all one has |i— Cjl ^ 
with equality only if Q e R. Hence 



\fiw)+tgiw)\ = \hiw)\ > \h{w)\ 



h{w) = \f{w) - igiw)\. 



To prove (3) we need to show that the univariate polynomial g(z„_|_i) = f{w) + 
Zn+ig{w) is stable. If g{w) = then q{zn^i) is a non-zero constant so it is stable. 
Therefore we may assume that g{w) ^ 0. By the above we have \ f{w)lg(w) + i| > 
\f(yS)lg(w) — i|, which implies that Im(/(u))/(7(u')) > 0. Now if q(2:„+i) — then 
f{w)/g{w) + Zn+i = and thus Im(z„+i) < 0. It follows that q G Hi(C) and since 
w S {C e C : Im(() > 0}" was arbitrarily fixed we deduce that f+Zn+ig € H„-(-i(C) 
hence / + Zn+ig € 7i„+i(R), which confirms (3). 

Finally, to show that (1) (2) note that by the above it is actually enough to 
prove that (3) (2). Now if (3) holds then letting Zn+i = i we get in particular 
that h £ Hn{C) and therefore \f{w) + ig{w)\ > \f{w) — ig{w)\ for all w G {C G C : 
Im(C) > 0}". It remains to show that we cannot have |/(u')+«g(?i') I — \ f{w)—ig{w)\ 
for such w. Supposing the contrary it follows from the above arguments that 
for some a G M" and (3 G M" all the zeros Cj: 1 ^ i ^ d, of the polynomial 
p{t) = /(" + tf3) + ig{a + tj3) are real. Since f{a + t/3) and g{a + tj3) have real 
coefficients this can occur only if these polynomials are constant multiples of each 
other, say /(a -I- tp) + ag{a + t(3) = for some a G R. By setting the variable 
Zn+i = a we see that f{z) + ag{z) is either a stable polynomial or identically zero. 
Since it vanishes for z — a -\- i(3 it must be identically zero, which contradicts the 
assumption that / and g are not constant multiples of each other. We conclude 
that (3) ^ (2), which completes the proof of the lemma. □ 

Remark 1.3. The equivalence (3) <^ (4) in Lemma [1.81 extends in obvious fashion 
to Hg-stable polynomials with complex coefficients: if /, g G C[zi, . . . , z„] \ {0} then 
/ + Zn+ig is Hg-stable (in n + 1 variables) if and only if 

Im('£!!M^>o, z=(zi,...,z„)GlHI^ 
V 9{z) J 
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Note in particular that for the open right-half plane the above relation reads 
Re( 44)^0, Re{zj) > 0, I < j <n. 



Definition 1.2. For g, f E R[zi, . . . , 2:„] and I < j < n define the j-th Wronskian 
of g, f as 

The following multivariate analog of the Hermite-Kakeya-Obreschkoff theorem 
(see, e.g., (43i Theorem 6.3.8] for the classical univariate version) was first proved 
in [S] . Below we give an alternative shorter proof. 

Theorem 1.9. Let f^gE . . • , Zn]. All non-zero polynomials in the space 

{af + f3g : a,/3 e K} 

are stable if and only if either f ~ g = 0, f g or g <ti f ■ Moreover, if g ^ f then 
Wj[g, f](x) < for all x G R" and 1 < j < n, and if f <. g then Wj[g, f](x) > 
for all X e R" and 1 < j < n. 

Proof. If / and g are linearly dependent then both conditions say that / and g are 
real stable or zero and then the Wronskians are all zero. Hence we may assume 
that /, 5 G R[zi, ■ ■ ■ , Zn] \ {0} are not constant multiples of each other. 

The "if" direction follows from Lemma [1.81 (1) ^ (3) and Lemma [T77] (1). 

For the converse assume that all non-zero polynomials in the space {af + f3g : 
a, /3 € R} are stable and that there exist z,w E {C E C : Im(C) > 0}" for which 

lJl^)>0 and lJl^)<0. 



By connectivity there is a number v E {( E C : Ini(C) > 0}" for which f{v)/g{v) — 
a G R \ {0}. It follows that / — ag is identically zero which we assumed was not 
the case. Hence Im(/(z)/g(z)) has constant sign for z G {C G C : Im(C) > 0}" and 
we conclude, by Lemma FOl (4) ^ (1), that either / <C g or g <C /. 

Suppose now that g f, where we may assume that g is not identically zero. 
Let X G K" be such that g{x) ^ and consider the rational function q{t) — 
f{x + ejt)/g{x + ejt), where Cj is the j-th standard basis vector of R". Hence q{t) 
is analytic at the origin with a first order Taylor expansion given by 



g{x) g^{x) 



From Lemma [LSI and the fact that lm{f{x)/ g{x)) = we then get Wj[g, f]{x) < 0. 
Since the set of all x E R" for which g{x) is dense in R" the theorem follows. □ 

Corollary 1.10. If f, g E ^[zi, . . . , z„] \ {0} are real stable polynomials such that 
f g and g <ti f then f — ag for some a G R. 

Proof, li f,g E ^[zi, . . . , z„] \ {0} are real stable polynomials such that f ^ g and 
g ^ f then f g and / ^ —g. Hence, by Theorem II. 9[ 



for all 1 < j < n and x E R" for which f(x) ^ 0. Thus g// is a real constant. □ 
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2. Hard-Core/ Algebraic Classification: Sufficiency 

2.1. Multi-afRne Polynomials and the Lieb-Sokal Lemma. A polynomial 
/ G C[zi, . . . , Zn] is multi-affine if all monomials in its Taylor expansion are square- 
free, i.e., if / can be written as 

f{z) = a{S)z^, where := [| z^, a{S) e C, [n] = {1, . . . , n}. 

sc[ii] ies 

Hence / G C[zi, . . . , z„] is multi-affine if and only if / G C(in)[zi, . . . , z„], where 
(1") = (1,...,1) gN". 

The following lemma (in the case n = 1) is due to Lieb and Sokal '33', Lemma 2.3] . 
For completeness we provide here a short proof essentially based on the same idea. 

Lemma 2.1 (Lieb-Sokal). Let P{z) + wQ{z) G C[zi, . . . , Zn,w] be stable. If the 
degree in the variable Zj is at most one then the polynomial 

is either identically zero or stable. 

Proof. We may assume that Q{z) is not identically equal to zero and that j = 1. 
Since Q{z) is stable and lm{w) > <;4> lm{—w^^) > the polynomial 

WQ[Z) ^■wQ(Zi-W ,Z2,...,Zn) 

ozi 



is stable. Hence, by Remark 11.31 (with 6 = 0) one has 

V ) \Qi^)J V J - 

for all z G {C G C : Im(C) > 0}" , so by Remark 1 1 . 31 again the polynomial 

Piz)^'-§^+wQiz) 

is stable. In light of Lemma 1 1.71 (1), this implies in particular the conclusion of the 
lemma (letting w = 0). □ 



We can now settle the sufficiency part of Theorem 1 1.1 1 in the case of multi-affine 
polynomials, i.e., for k = (1, . . . , 1). 

Lemma 2.2. Let T : C(ij.)[zi, . . . , z„] — )■ C[zi, . . . , z„] be a linear operator such 
that 

Gt(2,w) :=T[(z + u;)["lj = ^ T [z^] u.["l\^ 

5C[ri] 

is stable, where z = (zi, . . . , z^) and w = {wi, . . . , w„). Then T preserves stability. 
Proof. Let T be as in the lemma. Since lv[i{wj) > Im(— w^^) > we have that 

e ■H2n{'C) ^ Gt{z, w) g H2„(C). 
Therefore, if / G C[wi, . . . , u„] is stable and multi-affine we deduce that 
w^-^GT{z,-~w-')f{v) = ^ T [z^] {~wff{v) G H3„(C), 
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where v = {vi, . . . , By repeated use of Lieb-Sokal's Lemma 12711 we then get 
^ T[z''] /(^)(«)e7^2„(C)U{0}. 

SC[n] 

Letting {C G C : Im(C) > 0}" 3 v ^ and hivoking Hurwitz' Theorem 11.61 we 
obtain 

nf) - E ^ e u {0}, 

SC[n] 

which proves the lemma. □ 

2.2. Decoupling Schemes: Polarizations of Polynomials and Operators. 

For K — . . . , K„) G N" let be the space of multi-affine polynomials in the 

variables {z^ '■ ^ < i < n,! < j < Ki}. Define a (linear) polarization operator 

nT :C«[zi,...,z„] ^CJ,^ 

that associates to each / e Ck[zi, . . . , 2;„] the unique polynomial nj.(/) e C^,f^ 
such that 

(a) for any 1 < i < n the polynomial Il}.{f) is symmetric in {zij '■ 1 < j < Ki}; 

(b) if we let Zij = Zi for all 1 < i < n and 1 < j < Ki in n^(/) we recover /. 
In other words, if a < k then 

nK-Z") = Ea^izii, . . . , Zii^J ■ ■ ■ Ea,^{Zni, . . . , Zni^J, 

where Ei{xi, . . . , Xm) is the i-th elementary symmetric polynomial in the variables 
. . . , Xyji^ that is, 

Ei{xi,...,Xm)= E < I < m. 

SC[m], |S|=i 

Dually we define a (linear) projection operator 

ni:C^A^C«[zi,...,z„] 

by letting z^ i-^ and extending linearly. Note that by (b) above o Ii\ is the 
identity operator on [zi , . . . , z„] while 11^ o nj, is the operator on C^^^ that for 
each 1 < i < n symmetrizes all the variables in {zij : 1 < j < Ki} . 

Remark 2.1. In physics literature polarization and projection operators as above 
have mostly been used for univariate polynomials and are often referred to as "de- 
coupling procedures" [31 [ST] . Such procedures combined with Polya type results 
Hilfssatz II] were employed by Lee and Yang in their original proof of the "circle 
theorem" [3TJ Appendix II] (see also Kac's comments in pp. 424-426] and [SI 
§9]). Multivariate polarization notions similar to those defined above are used in 
the theory of hyperbolic polynomials and partial differential equations [HI US] . 

Let T : C k[zi, . . . , z„] — * C^[zi, . . . , Zrt\ be a linear operator. The polarization of 
T is defined as the linear operator 11 (T) : C^j^ '^'aia given by 

n(r) = nToToni, (2.1) 

and conversely we have 

r = n:^„n(r)onT. (2.2) 
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It is immediate from Lemma 11.71 (4) that the projection operator Ilj: preserves 
stabihty, but the remarkable fact is that so does Hj. . This is essentiahy the famous 
Gr ace- Walsh- Szego coincidence theorem I^OIISSIISH] albeit in a disguised form. For 
completeness let us state this theorem. 

Theorem 2.3 (Grace- Walsh-Szego). Let f be a symmetric multi-ajfine polynomial 
in n complex variables, let C he an open or closed circular domain, and let ^i, . . . 
be points in C. Suppose further that either the total degree of f equals n or that C 
is convex (or both). Then there exists at least one point S, € C such that 

= (2.3) 

We give a new and self-contained proof of Theorem 12.31 in [6, §2]. 

Proposition 2.4. Let f E Ck,[zi, . . . , Zn]- Then f is stable if and only ifll}.{f) is 
stable. 

Proof. The "if" direction follows from the fact that Ilj. preserves stability. For the 
other direction it suffices to prove the theorem for univariate polynomials since we 
may polarize the variables one at a time (when the other variables are fixed in the 
upper half-plane). Suppose therefore that f{z) = J2k=o ^^z'' is a univariate stable 
polynomial. Then the polarization of / is given by 

d / ,x -1 



F{zi,...,Zn) 



k=0 



Ekizi, 



,Zd) 



If F{C,i, . . . , C„) = for some (Ci, . . . , Cn) G then by the Grace- Walsh-Szego 
theorem (Theorem E^l) there exists C, E H such that /(() = F(C, . . . , C) = 0, 
contradicting the stability of /. □ 



Lemma 2.5. Let T 



zi , . . . , z„] be a linear operator. Then 



the symbol of the polarization of T is the polarization of the symbol of T, that is. 



Gn(T) - n^eftCCT), 



where Ujy^,., : C^^^zi, ■ 
Proof. We have that 



,Wl,. . . ,Wr, 



2=1 j=l 



Wij) 



(nI^onI-)[(z + ^)'=] 



where the upper indices z and w indicate that the corresponding operators act only 
on the z-variables and w- variables, respectively. Hence the symbol of the polarized 
operator may be written as 



Gn{T) 



(nT^oTon^-) 



(nT-oroni^)[(nI 



i=lj=l 

on! 



W) 



Using the fact that 11;^^ o 11^^ is the identity operator and that T and H^^ commute 
(since T only acts on the z-variables) we get 

Gum = (n^^ ° o t) [{z + wr] = n^g^ jgt), 
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as desired. □ 



Proof of Sufficiency in Theorem ] If Gt is stable then so is Gn(T) by Lemma [2.5l 
and Proposition l2.4l From Lemma [2T2l we deduce that n(T) preserves stabihty, and 
then by equation (|2.2p so does T since in view of Lemma 11.71 and Proposition 12.41 
both operators njj, and Hj. preserve stabihty. □ 

3. Algebraic Classification: Necessity in the Complex Case 

We will need the following simple lemma which makes it clear that there are 
plenty of stable polynomials. 

Lemma 3.1. Let f e C^[zi, . . . , Zn], where k = (ki,...,k„) e W\ and W = 
(Wi, . . . , Wn) (z {C (z C : Im(C) > 0}". Then for all e > sufficiently small one has 

{z + + e W„(C), where (z + = [zi + Wi)^' ■ ■ ■ (z„ + W„)«". 

Proof Set Y = (Im(T4^i), . . . , lm{Wn)) e W]^. For a < k we then have 



(z + wy 



< Y"-"" if z e {C e C : Im(C) > 0}" 



Expanding / in powers of z + T4^ we see that there exists eo > such that 

<- if ze{CeC:lm(C)>0}". 



\iz + W)-\ eo 

Hence {{z + W)" + ef{z)\ > \{z + W)'^\-e\f{z)\ > for aU z e {C € C : Im(C) > 0}" 
and e £ (0, eo), so that in particular (z + W)'^ + ef{z) S H„(C) for all such e. □ 

Remark 3.1. From the proof of Lemma l3.1l it follows that the topological dimension 
of the set Hn{C) n C„[zi, . . . , z„] (of all stable polynomials in n variables of degree 
at most k) equals (1 + ki) • • • (1 + k„). In particular, this explains the statement 
made in [11] §1] concerning the topological dimension of the set of so-called strongly 
Rayleigh probability measures on the Boolean algebra 2^"!, where [n] — {I, . . . ,n}. 

Lemma 3.2. Let V C K[zi, . . . , z„] be a K-linear space, where K = M or C 

(i) //K = M and every non-zero element of V is real stable then dimy < 2. 

(ii) //K = C and every non-zero element of V is stable then dimT^ < 1. 

Proof. We first deal with the real case. Suppose that there are three linearly in- 
dependent polynomials /i,/2 and /s in V. By the multivariate Hermite-Kakeya- 
Obreschkoff theorem (Theorem ll.Qp and the assumption on V these polynomials are 
mutually in proper position. Without loss of generality we may assume that /i <C /2 
and /i S> /a. Now consider the line segment in V given by £e = 6*73 + (1 — d)f2, 
< < 1, and note in particular that for any such 6 one has either /i <C £e or 
/i > £0. Set 7] = snp{d e [0, 1] : /i < £e}- Since /i < £0 and fi > £1 it follows 
from Hurwitz' theorem (Theorem 1 1.6p that /i ^ and /i ^ irj- This means that 
/i and £jj are constant multiples of each other (cf. CoroUarv ll.lOp . contrary to the 
assumption that /i, /2 and /a are linearly independent. 
For the complex case let 

Vr ^ {p : p + iq e V with p,q e ]R[zi, . . . , z„]} 

be the "real component" of V. By Theorem 11.91 all non-zero polynomials in Vr 
are real stable hence dimR Vr < 2 by part (i) proved above. Note that V is the 
complex span of Vr. If dimR Vr < 1 we are done so we may assume that {p, q} is 
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a basis of Vr with f := p + iq Cz V which by assumption is a (not identicaUy zero) 
stable polynomiaL By Theorem 11.91 again we have > 0, 1 < j < n, 

X G M", and since p and q are linearly independent the jo-th Wronskian g] 
is not identically zero for some index jq. Now if g £ V \ {0, f} we may write 

g — ap + bq + i{cp + dq) 

for some a, h,c,d £ M. We have to show that g is a (complex) constant multiple of 
/. Since g is not identically zero it is a stable polynomial, so that in particular 

Wjg [ap + bq, cp + dq] = {ad - bc)Wjg [p, q] > 

(cf. Theorem [LH). Since Wj^ [p, q] > and Wj^ [p, g] ^ it follows that ad-bc>0. 
Let M, u e R and note that by linearity we have 

g + {u + iv)f = (a + u)p + {b - v)q + i((c + v)p + (d + u)q) € V. 

Arguing as above we deduce that 

H{u, v) := (a + u){d + u) ~ {b - v){c + t;) > 

for aU u,v eR. But 

4:H{u, v) = (2u + a + d)^ + {2v + c- bf - (a - df - (6 + cf , 

so i?(u, w) > for all u, t; e R if and only if a = d and b = — c. This gives 

g ~ ap — cq + i{cp + aq) = (a + ic){p + iq) = (a + ic)/, 

as was to be shown. □ 

We can now show that the symbol of any stability preserver whose range is not 
one-dimensional must necessarily be stable. 

Proof of Necessity in Theorem ll.li Suppose that T : C„[zi, . . . , z„] — > C[zi, . . . , z„] 
preserves stability Given W & {( e C : Im(C) > 0}" we have that T[{z + W)""] G 
HniC) U {0}. Assume first that there exists W € {C € C : Im(C) > 0}" for which 
T[{z + W)"] EE and let / e C„[zi, . . . , z„]. Then by Lemma O there is an e > 
such that (z + W)" + ef{z) e UniC). It follows that 

eT[/(z)] = T[{z + WT + ef{z)] e Hn{C) U {0}. 

Hence the image of T is a linear space whose non-zero elements are all stable 
polynomials. By Lemma [321 (ii) we have that T is given by T{f) = a{f)P, where 
P is a stable polynomial and a is a linear functional. 

On the other hand, if T[{z + W)""] ^ for € {C € C : Im(C) > 0}" then 
T[{z + Wy] e n„{C) for all such M^, so that Gt(^,w) = T[iz + wy] € H2„(C). □ 

4. Real Stability Preservers 
The aim of this section is to prove Theorem 11.21 Let 

'H-iC)^{f{-z):f{z)enn{C)}. 
Proposition 4.1. For any n E N the following holds: 

HniC) n n-iC) = CW„(R) := {c/ : c e C,/ e W„(M)}. 
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Proof. Let h,g G M[zi, . . . , z„] be such that h + ig £ 'Wn(C) n (C) and set 
f = h-ig. lfx,ye R" then 

f{x + iy) ^ h{-{-x + iy)) + ig(-(-a; + iy)). (4.1) 

Since h + ig ^ 7Y,7(C), the right-hand side of (|4.1|) is non- vanishing whenever 
y S M" , which amounts to saying that f{x + iy) ^ for aU y S R" . By Definition 
ll.ll this imphes that —g <C /i, i.e., h ^ g. On the other hand, g h since h + ig G 
H,i(C). Therefore the proposition is just a reformulation of CoroUarv ll.lOl □ 

Proof of Theorem l 1 . Si If T : Rk[zi, . . . , 2;„] — > R[2;i, . . . , z„] is a hnear operator 
as in part (a) of Theorem 11.21 then it preserves real stability (by Theorem II. 9p . 
Moreover, if T, T' : R„[zi, . . . , z„] — > R[zi, . . . , Zn] are linear operators whose sym- 
bols satisfy Gt'{z, w) = Gt{z, —w) then these operators are related by T'{f){z) — 
(— l)T(/(— z)), hence they preserve real stability simultaneously. To settle the suf- 
ficiency part of the theorem we may therefore assume that Gt{z,w) is real stable. 
But then the desired conclusion simply follows from the complex case. 

To prove the necessity part of Theorem 11.21 consider Gt{z, w) — T[{z + w)'^] and 
letW £{zeC: Im(z) > 0}". Write (z + W)" G H„(C) as 

(z + wr - F{z) + iG{z), F,Ge R.[zi, . . . , z„]. 

By the multivariate Hermite-Kakeya-Obreschkoff theorem (Theorem 11.91) we have 
that aF + I3G e H„(R) U {0} for aU a,/3 € R and hence that aT{F) + (3T{G) e 
H„(R) U {0} for all a,/3 e R. By Theorem [T9l again this means that 

T{F) + iT{G) = Gt{z, W) e Hn{C) U H-{C) U {0}. 

Suppose that there are Wi,W2 e {z e C : Im(z) > 0}" such that Gt(z, W^i) e 
H„(C) U {0} and G't(z, W2) G H,T(C) U {0}. By a homotopy argument we then 
deduce that there exists t e [0, 1] such that 

GTiz,w') = HniC) nn-ic) u {0}, 

where W = {l-t)Wi+tW2 G {z G C : Im(z) > 0}". Therefore, by PropositiongT] 
there is a real stable polynomial p and o, 6 G R such that Gt{z, W) = {a + hi)p{z). 
Let us write (z + W')'^ — f{z) + ig{z), where f,g have real coefficients. Then 
T{f) ^ ap and T{g) = bp so that 

Tibf - ag) = 0. 

Now, as noted in the proof of Lemma 13.11 for any h G R^ [zi , . . . , z„] there is an 
e > such that 

\h\ 1 
\{z + W'Y\ e 

whenever Im(zi) > 0, 1 < z < n. Hence, if + 6^ 7^ (the case a = 6 = follows 
similarly) we get 

1^1 _ 1^1 / \(i ■ _ 1 

\bf-ag + i{af + hg)\ ~ \{b + ia){z + W'Y\ ^ ^''^ ^'"^'^ ~ ^' 

This implies that 

\bf -ag + e'h + i{af + bg)\>Q 
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whenever Ini(zi) > 0, 1 < i < In particular, bf — ag + e'h + i{af + bg) G 7i„(C) 
so by Theorem [L9l we have bf ~ ag + e'h e n„iR) U {0}. It foUows that 

T{h) = ^T{bf -ag + e'h) e 7^„(M) U {0}. 

This means that aU nonzero polynomials in the image of T are real stable and by 
Lemma 13.21 we conclude that the image of T is of dimension at most two. 

Thus we may assume that T[{z+W)''] G HniC) for aU G {z e C : Im(z) > 0}" 
or T[{z - W)'^] e n„iC) for edl W e {z e C : lm{z) > 0}". But this amounts to 
saying that r[(z + w)"] G H2„(K) or T[{z - w)"] G H2n(K), as claimed. □ 

5. Soft-Core/Transcendental Classifications 

In this section we settle Theorems 11.31 and 11.41 that is, the transcendental char- 
acterizations of stability, respectively real stability preservers. Let K = C or 
R. A linear operator T : IK[zi, . . . , z„] ]K[zi, . . . , z„] preserves stability if and 
only if the same holds for all its restrictions : ^k[zi, ■ ■ ■ , Zn] IK[zi, . . . , z„], 
K G N". Therefore, by Theorem 11.11 Theorem 11.21 and the definitions of the sym- 
bols Gt^ (z, w), Gt{z, w) we see that Theorem ll.31 respectively Theorcm ll.4l would 
follow from the case K = C, respectively IK = K, of the next theorem. As before 
we use the standard (product) partial order < on N": for a, (3 € N" one has a < f3 
if ai < Pi for all 1 < i < n. Set a! = ai! • • • a„! as in SjT]and let 

. ^ f ttAtt if a < /9, 
1 otherwise. 

Theorem 5.1. Let F{z,w) — X^aeN" where z — (zi,...,z„) and w — 

{wi, . . . , Wn), be a formal power series in w with coefficients in K[zi, . . . , z„]. Then 
F{z, w) G 7i2ri(K.) if and only if 

^(/3)„P„(z)u-" G W2n(IK)U{0} 
a</3 

for all /3 G N". 

The proof of Theorem [5Tl] requires several new ingredients and additional results 
that we proceed to describe. Since the arguments are the same for K = R and 
K = C, we will only focus on the latter case. 

5.1. Generalized Jensen Multipliers. For a,l3 £ N" let J{a,P) = (/3)q/3~" 
(we use the convention that fc^*^ = 1 for fc = 0). For n — 1 these are the so-called 
Jensen multipliers that are known to preserve (univariate) stability [14[ I32[ I43j . 
For each fixed /3 G N" the sequences {J{a, (3)}a<p and {{(3)a}a<p may be called 
multivariate (or generalized) Jensen multipliers in view of the following lemma. 

Lemma 5.2. Let (3 G N". The linear operators on C[zi, . . . , z„] defined by 

z" ^ J{a,(3)z°', aGN", 
z" ^ {f3)az", aeN", 



preserve stability. 
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Proof. Fix /9 e N" . Since the first operator is a composition of the second operator 
and a rescaling of the variables, it is enough to prove the lemma only for the second 
operator - call it T and denote by r„ its restriction to Ck[zi, ■ • • , Zn], where k G N". 
By Theorem 1 1.1 1 we need to show that 

Gt. (z, w) = T4iz + wr] - J2 Q 



is stable for all k However, 



J J \3 



1 



so the lemma amounts to showing that for any m, n £ N the univariate polynomial 



m - E 



3=0 



is real- rooted (then necessarily with all negative roots). To prove this note that 



where S is the operator on C[t] given by S* = (l + ^)"- To prove the well known 
fact that S preserves stability it suffices to consider the case n = 1. The symbol 
of S is S[{t + tu)™] = {m + t + w){t + w)™"-^, which is a stable polynomial. The 
desired conclusion now follows from Theorem IL21 □ 



5.2. Multivariate Szasz Principles. A remarkable fact about univariate stable 
polynomials, first noted by Szasz 55] and subsequently used by Edrei [18] . is that 
the growth of such a polynomial is controlled by its first few non-vanishing coeffi- 
cients. We extend Szasz' principle to higher dimensions and establish estimates for 
the growth of multivariate stable polynomials. 

Lemma 5.3 (Szasz). Suppose that f{z) = 1 + X]i=i '^i-^* — 11^=1 (1 + 0-^) stable. 
Then 

k 

ElOP<3KP + 2|a2|. 

j=i 

Proof. By assumption one has Im(^j) < 0, 1 < j < fe, hence 

k / k \ ^ 

EMe,)'< EMe,) =Im{a,r. 
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Since X]^=i ^| ~ ^ ^02 it follows that 



k k k 



= Re(a2_2a2) + 2Elm(0)' 

< Rc{al - 202) + 2Im(ai)2 < 3|aip + 2|a2| 



as claimed. 



□ 



Szasz used his lemma to obtain the following bound for the polynomial / in 
Lemma 15.31 

|/(z)| <exp(r|ai|+3r2|ai|2 + 3r2|a2|), \z\ < r. 



Remark 5.1. If the univariate polynomial / in Lemma [5.3l is assumed to have all real 
and negative zeros then '^f — ^i- Using the latter inequality in the following 

arguments one can get stronger bounds for the growth of multivariate real stable 
polynomials with all non-negative coefficients. 

As before, we let {ej}"^]^ be the standard basis in M". 

Lemma 5.4. Suppose that the polynomial 

/(z) = l+ E a{f3)zP eC[zi,...,Zn] 

\p\>o 

is stable and let 



1/2 



Then 



whe 



A= 3 El«(e')l +2E + 



Proof. Let us rewrite the above polynomial as f{z) — J2a ^(q^)-^" and suppose that 
/3 7^ is such that a{(3) ^ 0. Since the operator 1-^ J(a, a S N", preserves 

stability fLemma l5.2p the polynomial 

J(/3, P)a{(3)z^ + J{a, P)a{a)z" + 1 

n<\a\<\0\ 

is stable. Setting all variables equal to t we deduce that the univariate polynomial 
g{t) of degree k = \P\ given by 



^ 2<|Q|<fc 
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Ci H- 67" 



^ (l-lM)a(2e.) 



i ■■ Pi>2 



is stable. If we rewrite g{t) — X]i=o ~ 11^=1(1 + ^j^) then by Lemma [5.31 we 



get 



-, fe/2 



< 



which is the desired estimate. 

Theorem 5.5. Suppose that the polynomial 



□ 



|/3|>0 



is stable and let 



B = 2'" 



e - 1 



2.0210. 



*;i = l 



Then 



^il < 1 < « < n- !• < 5e<^'' , r > 0. 



Proof. We make use of the fohowing simple inequalities that may be proved by 
induction or from Stirling estimates. 



e"" < — < (671+ l)e-", n>0. (5.1) 
n" 

Let d{n,k) = E/3eN'M/3|=fc/3''//^! and note that |{/3 e N" : |^| = /c}| = for 
fc G N. By (|5.ip we have the following (rough) estimate: 

+ k — 1^ 



d{n,k) < 



Let A be as in Lemma 15.41 and apply the same lemma: 

max {|/(z)| : |z.| < ^, 1 < * < < E l/^l" 



J2d{n,k)k-''/^{Arf < 2"-^ ^ fc-'^-/2(2eAr)'=. 



/£ = 



k=0 



Set i? = 2eAr. By (|5.ip and the Cauchy-Schwarz inequality we get 



J2k''^/'R'^<J2^{ek + l) 



k=0 



k=0 



IT 



< 



.^(efc + l)f 
\ \fc=o 



(5.2) 



\k=0 



1 . pflV2 



e- 1 



which combined with (|5.2p proves the theorem. 



□ 
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Given f{z) = Y.a a(a)2" G . . . , z„] let supp(/) = {a e N" : a{a) ^ 0} be 
its support. Denote by Ai{f) the set of minimal elements of supp(/) with respect 
to the partial order < on N". Moreover, for a fixed set A4 C N" we let 

= {a + /3 : a e G N", |/3| < £}. 

We can now establish a multivariate Szasz principle. 

Theorem 5.6. Let A4 C N" be a finite non-empty set and f{z) — X]q ^(o^)-^" ^ 
C[2;i, . . . , z„] be a stable polynomial with M{f) = M. Then there are constants B 
and C depending only on the coejficients a{a) with a Cz M2 such that 



max||/(z)| : \z,\ < r, 1 < i < n| < Be'^''' 



for all r > 0. Moreover, B and C can be chosen so that they depend continuously 
on the aforementioned set of coefficients. 

Proof. We construct B and C inductively over m = n + fc, where 

k = max{|a| : a £ M{f)}. 

If /(O) 7^ then fc = and the bound follows from Theorem [531 so we may assume 
that A; > 0. Also, if n = 1 we may just factor out a multiple of z and use the fact 
that r < e*" together with Theorem 15.51 Hence we may assume that n > 1. 

Let a — {ai, . . . , q;„) € A^(/) be such that \a\ = k and suppose (without loss of 
generality) that ai > 0. Now 

r-l Qf 



SO that 



/ df 

/(Zl, . . . ,Z„) = /(0,Z2, . . . ,Z„) + Zl J — (Zlt, Z2, 

df 



, Zn)dt 



|/(zi, . . .,Zn)\r < 1/(0, Z2, • ■ ■ , Zn)\r 



|/(0,Z2,...,z„)|r <-Bie^i'^ and 



< B.e 



C2r^ 



dzi 

where \g\r = max ||5(z)| : |zi| < r, 1 < i < n| for g e C[zi, . . . , Zn]. Since 
max{|/3| : (3 G M{df/dzi)} 

and /(O, Z2, • ■ • , z„) is either identically zero or a stable polynomial in n— 1 variables, 
by the induction assumption we have 

K 

dzi 

Set C = max{Ci, C2 + 1} and B = 2max{_Bi,i?2}- From the above estimates we 
get \ f\r < Be'^'^ , which completes the induction step. □ 

5.3. Proof of Theorem 15.11 We can now settle Theorem 15.11 and thereby - as 
explained at the beginning of this section - Theorems 11.31 and 11.41 as well. Note 
that the case K = K in Theorem 15.11 follows from the case K = C since 7^271 (C) n 
R[zi, . . . , z„, wi, . . . , Wn] = H„(M). Therefore we assume that IK = C. 

Suppose that F{z,w) = SqgN" ^a(^)^" ^ W2«(C) and P„ e C[zi, . . . , z„], 
a G N". If Fjn{z,w) = X^qgN" Prn.a{z)w°' IS a scqucnce of polynomials in H2n(C) 
converging to F(z,w), uniformly on compacts, then Pm,a{,z) Pa{z) a,s m ^ 
00 uniformly on compacts for fixed a G N". Given (3 G N" let A be the linear 
operator on C[wi, . . . ,Wn] defined by 1-^ {P)aW°', a G N", and extend it to 
a linear operator A on <C[zi, . . . , Zn,wi, . . . ,Wn] by letting A[z'^w°'] = z'^A[w°'], 
a,-f £ N". We claim that A preserves stability in 2n variables. Indeed, this amounts 
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to saying that its restriction to Ck[zi, . . . , z„, wi, . . . , Wn] preserves stability for 
each K e N^". The latter statement is an immediate consequence of the fact that 
A preserves stability in n variables (by Lemma 15. 2p combined with Theorem 11.11 
and a straightforward computation showing that the symbol of A^ is a stable 
polynomial in An variables. Now, since A sends all but finitely many monomials w"' 
to we have that A[Fm{z,w)] A[F{z,w)] as m — > oo, uniformly on compacts, 
which gives A[F{z, w)] S H2n{^) U {0}. This proves the necessity part of Theorem 

To prove the converse assume that J2a<f3if^)aPa{z)w°' € H2n(C) U {0} for all 
/3 G N". Let TO e N, set (3m — {m, . . . , m) and note that 

Fm{z,w) J2 J{a,P.n)Pc.{z)w'' = iPrn)Mz) (^)" G H2„(C) U {0}. 

We want to show that there is a positive integer M such that for all r > there is a 
constant Kr (depending only on r) such that for all to > M we have \Fm {z,w)\ < Kr 
whenever \wj\ < r and \zj\ < r for 1 < j < n. This will prove the sufficiency part 
of Theorem 15.11 since {Fm{z,w)}meN is then a normal family whose convergent 
subsequences converge to F{z,w) (by the fact that limm^oo >^(aj An) — 1 for all 
a G N"). Now, there exists a positive integer N such that the set A^(F,„) of 
minimal elements of supp(F„i) does not change for m > N (this is easily seen from 
the fact that all lower degree terms are multiplied by eventually non-zero numbers 
J{a, (3m))- Since limm^oc J {ct, (3m) = 1 it follows from Theorem 15.61 that there is 
a positive integer M and constants B and C such that if \wj\ < r and \zj\ < r, 
1 < J < n, then 

\Fmiz,w)\< Be'^''\ m>M, 
which completes the proof of the theorem. 

6. Circular Domains and Lee- Yang Polynomials 

We will now solve Problems [T] [2] whenever is a product of arbitrary open circu- 
lar domains. Using this we then characterize all linear transformations preserving 
the class of Lee- Yang polynomials defined with respect to any such set VI. 

6.1. Products of Open Circular Domains. Recall that a Mobius transforma- 
tion is a bijective conformal map of the extended complex plane C given by 

4>{C,) = a, 6, c, c? € C, ad — 6c = 1. (6.1) 

Note that one usually has the weaker requirement ad — bc^O but this is equivalent 
to (j6.ip which proves to be more convenient. An open circular domain is the image 
(in C) of the upper half-plane H under a Mobius transformation, i.e., an open disk, 
the complement of a closed disk or an open affine half-plane. 

As usual, the support of a polynomial f{z) = J2aeti" o-{ce)z°' g C[zi, . . . , z„], i.e., 
the set {a G N" : a{a) ^ 0}, is denoted by supp(/). We say that / is of degree 
K — (ki, . . . , K„) £ N" if deg^. (/) — Hi, i S [n], where as before [n] = {1, . . . , n}. 

To deal with discs and exterior of discs we need some auxiliary results. 

Lemma 6.1. Let {Ci}f^i be a family of circular domains, f e C[zi, . . . , z„] be of 

degree k G N", and J C [n] a (possibly empty) set such that Cj is the exterior of 
a disk whenever j e J. Denote by g the polynomial in the variables Zj, j G J, 
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obtained by setting Zi ~ a G Ci arbitrarily for i ^ J. If f is Ci x ■ ■ ■ x Cn-stable 
then supp(g) has a unique maximal element 7 with respect to the standard partial 
order on N'^ . Moreover, 7 is the same for all choices of Ci € Ci, i ^ J. 

Proof. Let us first prove the lemma in the case when J = [n]. We may assume 
that Cj is the exterior of the closed unit disk D for all j S [n] . Suppose that there 
is no unique maximal element of supp(/). For fixed A = (Ai, . . . , A„) G R!f. and 
(3 — (/3i, ...,/?„)€ N" consider the univariate polynomial 

F{X,/3;t) ■.^Kf{X,t''\...,Xy-), 

where K — K{X, /3) G C is the appropriate constant making F monic. Let A{(3) be 
the closure (with respect to uniform convergence on compacts) of the set {F{X, (3; t) : 
X G (1,00)"}. By Hurwitz' theorem (Theorem [m]) all elements of A{f3) are C \ D- 
stable. However, since there is no unique maximal element of supp(/) there exists 
P G N" such that A{f3) contains polynomials of different degrees. Hence, as we vary 
A at least one zero must escape to infinity. This is impossible since all the zeros of 
F{X,/3;t) have modulus at most one. 

If J ^ [n] and the maximal elements are different for different choices of Ci then 
as above some zero of a specialization of / would escape to infinity contradicting 
the boundedness of the zeros. □ 

When dealing with non-convex circular domains one has to be careful with de- 
grees. We are therefore led to consider the following classes of polynomials. 

Notation 6.1. If {Ci}^^^ is a family of circular domains and k = (ki, . . . , k„) G 
N" we let A/'k(Ci, . . . , C„) denote the set of Ci x • • • x C„-stable polynomials in 
Ck [^1 , . . • , Zn] that have degree kj in zj whenever Cj is non-convex. (If all Cj are 
convex then Af^(Ci, . . . ,Cn) is just the set of all Ci x • • • x C„-stable polynomials 
in Ck[zi, . . . , Zn] ) Let also 



AA(Ci,...,C„) = U AA,(Ci,...,C„). 



Lemma 6.2. Suppose that Ci, . . . , C„, Z3i, . . . , -D„ are open circular domains and 
K = (ki, . . . , Kn) G N". Then there are Mobius transformations 

as in (|6.ip such that the (invertible) linear transformation ; Ck[zi, . . . , z„] — > 
<Ck.[z\,..., Zn] defined by 

$„(/)(zi, . . . , Z„) = (ciZi + diT' ■ ■ ■ {CnZn + dn)^"f{Mzi), M^n)) (6.3) 

restricts to a bijection between ^^{Ci, . . . , C„) and JV^iDi, . . . , Z3„). 

Proof. The relation between tuples of open circular domains defined by relating 
tuples for which there is such a linear transformation <I>k is an equivalence relation. 
Therefore it suffices to assume that Cj = Dj for j > 2 and that Di = H {= Mq). 
Thus we want to find $k of the form 

$„(/)-(cZi+rf)-vf^^,^2,...,. 

V czi + d 
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If Ci is convex we may chose <j){z) = {az + b)/{cz + d) as in (16. 1|) so that 4>{H) = Ci 
and cz + d 7^ for aU z (z H. If Ci is non-convex we may chose (j) so that 
(j){H \ {rj}) — Ci for some specific point rj H , namely rj = —d/c. 

Now let / e AA„(Ci, . . . ,C„) and e = (6, • • ■ , Cn) e ^1 x • • ■ x A.- If +d ^ 
then 4>{£,i) £ Ci and <I'k(/)(^) 7^ by construction. If c^i + d = then Ci is non- 
convex so by Lemma l6. II the polynomial zi 1— > f{zl,S,2^(,s,^ ■ ■ ■ iCn) may be written 
as X^fcLo o^fe^i with a„j 7^ 0. But then 



where (f){z) — {az + b)/{cz + d) is chosen as above. It follows that — cz + d 7^ 
for all z G Ci so it only remains to prove that the degree of $^^(5) in zi is ki 
when Ci is non-convex. To do this fix £ Dj, 2 < j < n, and consider the 
_ff-stable polynomial g(zi, ^2, • ■ ■ , ^n) — Sfc=o "^fe-^i ■ The coefficient in front of z^^ 
in (f>^^((7)(zi, ^2, ■ • • , Cn) is (— c)"! X^fclo which is non-zero since as we 
previously noted one has —d/c E H. □ 

The following theorem answers a more precise version of Problem [T] for K = C 
and Q = Ci X ■ ■ ■ X C„, where the C; are arbitrary open circular domains. 

Theorem 6.3. Let k e N" and T : Ck[zi, . . . , Zn] C[zi, . . . , z„] he a linear op- 
erator. Let further d, i € [n], be open circular domains given by d — 
i £ [n], where 4>i, i G [n], are Mobius transformations as in (|6.2p such that the cor- 
responding linear transformation $k defined in (|6.3p restricts to a bijection between 
N^{H,..., H) and 7V;(Ci, . . . , C„) (cf. Lemma\6^. Then 



where a is a linear functional on [zi , . . . , z„] and P is a Ci x ■ ■ ■ x Cn - 
stable polynomial, or 
(b) the polynomial in 2n variables Zi, . . . , z„, Wi, . . . , Wn given by 



.4=1 

is Ci X ■ ■ ■ X Cn X Ci X ■ ■ ■ X Cn-stable. 

Remark 6.1. By analogy with the definitions made in m.l\ we may call the poly- 
nomial in Theorem 16. 31 fb) the algebraic symbol ofT with respect to Ci x ■ ■ ■ x Cn- 
Note that in the case of the open unit disk Ci — D, i £ [n], or the exterior of the 
closed unit disk Ci = C \ D, i e [n], this symbol becomes a constant multiple of 
T[(l -\- zw)'^] while in the case of an open half-plane bordering on the origin it is 
just a constant multiple of Gt{z,w) = T[(z -I- wY]. For the open right half-plane 
Mil it is often more convenient (but equivalent) to choose the symbol T[(l + zwY]. 



Hence e A^^I^i, . . . , 

If (/ e Nk,{Di, . . . , Dn) consider 



a^^{ad - bcY^ I c""^ ^ 0. 





n 



T {{aiZ^ + bi){c^w^ + di) + {aiWi + bi){ciZ^ + d^)) 
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Proof of Theorem \6.3i Let T : Ck,[zi, . . . , Zn] — > C[zi, . . . , Zn] be a linear operator 
and 7 = (71,..., 7„) e N" be such that T{Cn[zi, . . . , Zn]) C C-y[zi, . . . , Zn]. The 
theorem is trivial if the range of T has dimension at most one so we may assume 
that this is not the case. With (f>i, i e [n], and as in the statement of the theorem 
let : <Cj [zi, . . . , Zn] — > C-y [zi, . . . , Zn] be the linear transformation defined as in 
(fO]) only with 7 = (71, . . . , 7„) instead of k = . . . , k„). Set T' := $~iT$„. 
By Theorem 11.11 and Lemma 16.21 we have the following chain of equivalences: 

T(AA«(Ci, . . . , C„)) c AA(Ci, . . . , C„) U {0} ^ 

T' preserves (i7-)stability <;=^> 

G(z,w) := r'[(z-|-w)''] is (i/-)stable ^ 

$^0k(G) is Ci X • • • X C„ X Ci X • • • X C„-stable, 

where $^0^ : C^e^izi, . . . , z„, wi, . . . , w„] ^ C^e„[zi, . . . ,z„,wi, . . . ,w„] is given 

by 

■■■,Zn,Wi,.. .,Wn) 

n n 

= ^(CiZi + di^' Y\_i<^i'^i + di)'''fi(f>lizi), ■ ■ . ,(f>n{Zn),(f>l{wi), . . . ,0„(w„)), 
i=l 1=1 

see ()6.3p . Now $^0k(G) is precisely the polynomial in Theorem 16.31 (b). Indeed, 
letting subscript z and subscript w indicate that the corresponding operator acts 
on the z and w variables, respectively, we have 



which completes the proof of the theorem. □ 



Note that if all Ci are open convex circular domains then Theorem 16.31 answers 
Problem [T] (for such domains) precisely as stated in the introduction, and when 
non-convex circular domains are involved it answers a more precise version of the 
problem. 

6.2. Linear Preservers of Lee- Yang Polynomials. Given an open circular do- 
main G C C let G"" be the interior of the complement of G, that is, G*" = C \ G. 

Definition 6.1. Let k = . . . , k„) e N" and {G;}"^]^ be a family of open circular 
domains. The set of k- Lee- Yang polynomials with respect to {Gi}"^]^ is 

cy^iCi, . . . , G„) - AA,(Gi, . . . , G„) n AA,(G[, . . . , g;). 

A Lee- Yang polynomial with respect to {Ci}^^i is an element of the set 
£y{Ci,...,Cn):^ U £3^.(Gi,...,G„). 

This definition is motivated by the original Lee- Yang theorem [31] and its various 
versions such as the so-called "circle theorem" ([6l Theorem 8.4], [23], [49]). For 
instance, the class of "complex Lee- Yang polynomials" related to the latter which 
was introduced and studied in [3] corresponds to the case when Gi = D, i G [rt], 
in Definition 16. II Note also that by Corollarv 1 1 . 101 the Lee- Yang polynomials with 
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respect to where Hi — H , i E [n], are precisely the polynomials that are 

non-zero complex constant multiples of real stable polynomials. 

We say that a linear operator T : Ck[zi, . . . , z„] ^ C[zi, . . . ,Zn] preserves the k- 
Lee-Yang property with respect to {Ci}2^i if T{f) € Ly{Ci, . . . , C„)U{0} whenever 
/ G Cy i^{Ci, . . . , C„). From our symbol theorems we can deduce a characterization 
of all such linear operators. We only display the case when T is non- degenerate in 
the sense that its image has dimension larger than 2. 

Theorem 6.4. Let k e N", T : Ck[zi, . . . , Zn] — > C[zi, . . . , Zn] be a non-degenerate 
linear operator and Ci, i G [n], he open circular domains given by Ci = 4>'^^{H), 
i & [n], where (j^i, i € [n], are Mobius transformations as in Theorem \6.!A Then T 
preserves the n-Lee-Yang property with respect to {Ci\^^i if and only if either 

(a) the polynomial in 2n variables zi, . . . , Zn,wi, . . . ,Wn given by 



T 



{{a.A + hi){ciWi + di) + {a^Wi + bi){ciZi + d^))"' 



.i = l 



is Ci X • • • X C„ X Ci X • • • X C'n-stable and C[ x • • • x x C[ x • • • x C^-stable, 



(b) the polynomial in 2n variables zi, . . . , z„, wi, . . . , Wn given by 



T 



((a^Zj + bi){ctWi + d^) - {oiWi + bi){ciZ^ + d^)Y" 



i=l 



is Ci X • • • x C„ X Ci X • • • X Cn-stable and C[ x • • • x C,''j x C[ x • • • x C^^-stable. 

Proof. The proof of the corresponding theorem in the case n = 1 from [7] , namely 
[71 Theorem 3], extends naturally to the above setting. □ 



7. Related Results 

A multivariate polynomial is called strictly stable if it is non- vanishing whenever 
all variables lie in = {z G C : Im(z) > 0}. The arguments in fj2] carry over to 
the closed upper half-plane H and yield the following sufficient condition for strict 
stability preserving. 

Theorem 7.1. Let k G N" and T : Ck[zi, . . . , z„] — > C[zi,...,z„] he a linear 
operator. If Gt{z, w) is a strictly stable polynomial then T preserves strict stability. 

One can also check that the sufficiency part of Theorem 16.31 carries over to 
closed circular domains and thus yields a (sufficient) condition for linear operators 
T : Ck[zi, . . . , Zn] — ^ C[zi, . . . , z„] to preserve Ci x • • • x C„-stability, where Ci, 
i G [n], are arbitrary (open) circular domains. For simplicity we only state here a 
special case that extends Theorems 17. II 

Theorem 7.2. Let k G N", T : C„[zi, . . . , z„] C[zi,...,z„] a linear operator 
and C a convex closed circular domain given by C = (j)~^{H), where (j) is a Mobius 
transformation as in ()6.ip . // 

T [((az -t- h){cw + d) + {aw + b){cz + d)Y] 

is C-stable then T preserves C-stability. 
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As noted in Remark 16.11 for the the closed unit disk C = D the symbol used in 
Theorem l7.2l is just a constant multiple of T[{1 + zw)'^] while for a closed half-plane 
bordering on the origin it is a constant multiple of Gt{z,w) = T[{z + w)'^]. Note 
though that the conditions in Theorems 17. 1[ 17.21 are not necessary. For instance, 
the identity operator obviously preserves strict stability but its (algebraic) symbol 
{z + w)'^ is not strictly stable, cf. (i) in fJS] 

8. Further Directions 

To conclude we mention some of the most appealing cases where Problems [T]-[2] 
remain open: 

(i) is a closed circular domain, to begin with = H (cf. Theorem 17. ip : 

(ii) is a strip, e.g., = {z e C : |Im(2)| < a}, a > 0; 

(iii) is a sector, e.g., ft — {z^C: ^ — a< arg(z) < ^ + a}, a e (O, f ); 

(iv) is a closed or open strip or (double) sector. 

Finally, from a complex analytic perspective and its many potential applications 
[TU] it would be quite useful to establish analogs of these results for transcendental 
entire functions in one or several variables. 

Acknowledgments 

We wish to thank the Isaac Newton Institute for Mathematical Sciences, Uni- 
versity of Cambridge, for generous support during the programme "Combinatorics 
and Statistical Mechanics" (January- June 2008), where this work, [B] and [TT] were 
presented in a series of lectures. Special thanks are due to the organizers and the 
participants in this programme for numerous stimulating discussions. We would also 
like to thank Alex Scott, Bela BoUobas and Richard Kenyon for the opportunities to 
announce these results in June- July 2008 through talks given at the Mathematical 
Institute, University of Oxford, the Centre for Mathematical Sciences, University of 
Cambridge, and the programme "Recent Progress in Two-Dimensional Statistical 
Mechanics" , Banff International Research Station, respectively. 

References 

[1] T. Asano, Theorems on the partition functions of the Heisenberg ferromagnets, J. Phys. Soc. 

Japan 29 (1970), 350-359. 
[2] M. F. Atiyah, R. Bott, L. Garding, Lacunas for hyperbolic differential operators with constant 

coefficients I, Acta Math. 124 (1970), 109-189. 
[3] B. Beauzamy, On complex Lee and Yang polynomials, Comm. Math. Phys. 182 (1996), 177— 

184. 

[4] M. Biskup, C. Borgs, J. T. Chayes, R. Kotecky, Partition function zeros at first-order phase 
transitions: Pirogov-Sinai theory, J. Statist. Phys. 116 (2004), 97—155. 

[5] M. Biskup, C. Borgs, J. T. Chayes, L. J. Kleinwaks, R. Kotecky, Partition function zeros at 
first-order phase transitions: A general analysis. Comm. Math. Phys. 251 (2004), 79-131. 

[6] J. Borcea, P. Branden, The Lee-Yang and Polya-Schur programs. IL Theory of stable poly- 
nomials and applications, arXiv:0809.3087 

[7] J. Borcea, P. Branden, Polya-Schur master theorems for circular domains and their bound- 
aries, Ann. of Math., to appear; arXiv:math/0607416 

[8] J. Borcea, P. Branden, Applications of stable polynomials to mixed determinants: Johnson's 
conjectures, unimodality, and symmetrized Fischer products, Duke Math. J. 143 (2008), 
205-223. 

[9] J. Borcea, P. Branden, Multivariate Polya-Schur classification problems in the Weyl algebra, 
arXiv:math/0606360 



26 J. BORCEA AND P. BRANDEN 

[10] J. Borcea, P. Branden, G. Csordas, V. Vinnikov, Polya-Schur-Lax problems: hyperbolicity 
and stability preservers, Workshop Report, American Institute of Mathematics, Palo Alto, 
CA, May— June 2007; "http:/ /www. aimath.org/pastworkshops/polyaschurlax. html 

[11] J. Borcea, P. Branden, T. M. Liggett, Negative dependence and the geometry of polynomials, 
J. Amer. Math. Soc. 22 (2009), 521-567; arXiv: 0707. 2340 

[12] P. Branden, Polynomials with the half-plane property and matroid theory. Adv. Math. 216 
(2007), 302-320. 

[13] Y. Choe, J. Oxley, A. D. Sokal, D. G. Wagner, Homogeneous multivariate polynomials with 
the half-plane property, Adv. Appl. Math. 32 (2004), 88-187. 

[14] T. Craven, G. Csordas, Composition theorems, multiplier sequences and complex zero de- 
creasing sequences, in "Value Distribution Theory and Its Related Topics", G. Barsegian, I. 
Laine, C. C. Yang (Eds.), pp. 131-166, Kluver Press, 2004. 

[15] T. Craven, G. Csordas, W. Smith, The zeros of derivatives of entire functions and the Polya- 
Wiman conjecture, Ann. of Math. (2) 125 (1987), 405-431. 

[16] T. Craven, G. Csordas, Multiplier sequences for fields, Illinois J. Math. 21 (1977), 801-817. 

[17] G. Csordas, Linear operators and the distribution of zeros of entire functions. Complex Var. 
Elliptic Equ. 51 (2006), 625-632. 

[18] A. Edrci, Power series having partial sums with zeros in a half-plane, Proc. Amer. Math. 
Soc. 9 (1958), 320-324. 

[19] S. Fisk, Polynomials, roots, and interlacing. Versions 1-2, http://www.bowdoin.edu/fisk/, 
xx+700 pp. 

[20] J. H. Grace, The zeros of a polynomial, Proc. Cambridge Philos. Soc. 11 (1902), 352-357. 
[21] L. Garding, An inequality for hyperbolic polynomials, J. Math. Mech. 8 (1959), 957—965. 
[22] O. J. Heilmann, E. H. Lieb, Theory of monomer-dimer systems. Comm. Math. Phys. 25 
(1972), 190-232. 

[23] A. Hinkkanen, Schur products of certain polynomials, in "Lipa's Legacy: Proceedings of the 

Bers Colloquium", J. Dodziuk, L. Keenin (Eds.), Contemp. Math. 211, Amer. Math. Soc, 

Providence, RI, 1997, pp. 285-295. 
[24] L. Hormander, On a theorem of Grace, Math. Scand. 2 (1954), 55-64. 
[25] L. Hormander, Notions of Convexity. Progr. Math. 127. Birkhauser, Boston, MA, 1994. 
[26] A. Iserles, S. P. N0rsett, E. B. Saff, On transformations and zeros of polynomials. Rocky 

Mountain J. Math. 21 (1991), 331-357. 
[27] R. Kcnyon, A. Okounkov, S. Sheffield, Dimers and amoebae, Ann. of Math. (2) 163 (2006), 

1019-1056. 

[28] E. Laguerre, Fonctions du genre zero et du genre un, C. R. Acad. Sci. Paris 95 (1882), 
828-831. 

[29] P. D. Lax, Differential equations, difference equations and matrix theory. Comm. Pure Appl. 

Math. 6 (1958), 175-194. 
[30] G. F. Lawler, O. Schramm, W. Werner, Values of Brownian intersection exponents. I. Half- 
plane exponents. Acta Math. 187 (2001), 237-273; //. Plane exponents, ibid., 275-308. 
[31] T. D. Lee, C. N. Yang, Statistical theory of equations of state and phase transitions. H. 

Lattice gas and Ising model, Phys. Rev. 87 (1952), 410-419. 
[32] B. Ja. Levin, Distribution of zeros of entire functions. Transl. Math. Monogr. Vol. 5, Amer. 

Math. Soc, Providence, R.I., 1980. 
[33] E. H. Lieb, A. D. Sokal, A General Lee-Yang Theorem for One- Component and Multicom- 

ponent Ferromagnets, Comm. Math. Phys. 80 (1981), 153-179. 
[34] T. M. Liggett, Distributional Limits for the Symmetric Exclusion Process, to appear in Stoch. 

Proc. Appl., arXiv:0710.3606 
[35] M. Marden, The Geometry of the Zeros of a Polynomial in a Complex Variable. Math. Surveys 

Vol. 3, Amer. Math. Soc, New York, NY, 1949. 
[36] C. M. Newman, Gaussian correlation inequalities for ferromagnets, Z. Wahrscheinlichkeits- 

theorie und Verw. Gebiete 33 (1975), 75-93. 
[37] C. M. Newman, Inequalities for Ising models and field theories which obey the Lee-Yang 

theorem. Comm. Math. Phys. 41 (1975), 1-9. 
[38] C. M. Newman, Zeros of the partition function for generalized Ising systems. Comm. Pure 

Appl. Math. 27 (1974), 143-159. 
[39] G. Polya, Bemerkung iiber die Integraldarstellung der Riemannsche ^-Funktion, Acta Math. 

48 (1926), 305-317. 



THE LEE- YANG AND POLYA-SCHUR PROGRAMS. I. 



27 



[40] G. Polya, Collected Papers, Vol. II: Location of Zeros. Mathematicians of our Time. Vol. 8 

(R. P. Boas, ed.). The MIT Press, Cambridge, MA, 1974. 
[41] G. Polya, I. Schur, Uber zwei Arten von Faktorenfolgen in der Theorie der algebraischen 

Gleichungen, J. Reine Angew. Math. 144 (1914), 89-113. 
[42] G. Polya, G. Szego, Problems and Theorems in Analysis. Vol. II, Springer, Berlin, 1976. 
[43] Q. I. Rahman, G. Schmeisser, Analytic Theory of Polynomials. London Math. Soc. Monogr. 

(N. S.) Vol. 26, Oxford Univ. Press, New York, NY, 2002. 
[44] D. Ruelle, Zeros of graph- counting polynomials, Comm. Math. Phys. 200 (1999), 43-56. 
[45] D. Ruelle, Grace-like polynomials, Foundations of computational mathematics (Hong Kong, 

2000), pp. 405-421, World Sci. Publ. Co. Inc., River Edge, NJ, 2002. 
[46] D. Ruelle, Is our mathematics natural? The case of equilibrium statistical mechanics, Bull. 

Amer. Math. Soc. (N.S.) 19 (1988), 259-268. 
[47] D. Ruelle, Some remarks on the location of zeroes of the partition function for lattice systems. 

Comm. Math. Phys. 31 (1973), 265-277. 
[48] D. Ruelle, Extension of the Lee-Yang circle theorem, Phys. Rev. Lett. 26 (1971), 303—304. 
[49] D. Ruelle, Statistical mechanics: Rigorous results. Reprint of the 1989 edition. World Sci. 

Publ. Co. Inc., River Edge, NJ, 1999. 
[50] I. Schur, Zwei Sdtze iiber algebraische Gleichungen mit lauter reellen Wurzeln, J. Reine 

Angew. Math. 144 (1923), 75-88. 
[51] A. D. Scott, A. D. Sokal, The repulsive lattice gas, the independent- set polynomial, and the 

Lovdsz local lemma, J. Stat. Phys. 118 (2005), 1151-1261; arXiv: cond-mat / 0309352 
[52] S. Smirnov, W. Werner, Critical exponents for two-dimensional percolation. Math. Res. Lett. 

8 (2001), 729-744. 

[53] A. D. Sokal, The multivariate Tutte polynomial (alias Potts model) for graphs and matroids, 
in "Surveys in Combinatorics, 2005" (B. S. Webb, ed.), Cambridge Univ. Press, Cambridge, 
2005; arXiv:math.CO/0503607 

[54] A. D. Sokal, Chromatic roots are dense in the whole complex plane, Combin. Probab. Comput. 
13 (2004), 221-261. 

[55] O. Szasz, On sequences of polynomials and the distribution of their zeros. Bull. Amer. Math. 

Soc. 49 (1943), 377-383. 
[56] G. Szego, Bemerkungen zu einem Satz von J. H. Grace iiber die Wurzeln algebraischer 

Gleichungen, Math. Z. 13 (1922), 28-55. 
[57] D. G. Wagner, Weighted enumeration of spanning subgraphs with degree constraints, preprint 

arXiv:0803.1659 

[58] J. L. Walsh, On the location of the roots of certain types of polynomials. Trans. Amer. Math. 

Soc. 24 (1922), 163-180. 
[59] C. N. Yang, T. D. Lee, Statistical theory of equations of state and phase transitions. I. Theory 

of Condensation, Phys. Rev. 87 (1952), 404-409. 

Department of Mathematics, Stockholm University, SE-106 91 Stockholm, Sweden 
E-mail address: juliusOmath. su. se 

Department of Mathematics, Royal Institute of Technology, SE-100 44 Stockholm, 
Sweden 

E-mail address: pbrandenOmath. kth. se 



